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Abstract 

We develop a formalism to study tetraquarks using the generalized flip-flop potential, 
which include the tetraquark potential component. Technically this is a difficult problem, 
needing the solution of the Schrodinger equation in a multidimensional space. Since the 
tetraquark may at any time escape to a pair of mesons, here we study a simplified two- 
variable toy model and explore the analogy with a cherry in a glass, but a broken one 
where the cherry may escape from. We also compute the decay width in this two-variable 
picture, solving the Schrodinger equation for the outgoing spherical wave. 



1 Introduction, tetraquarks with flux tubes 



Our main motivation is to contribute to understand whether exotic hadrons exit or not. 
Although there is no QCD theorem ruling out exotics, they are so hard to find, that 
many friends even state that either exotics dont exist, or that at least they should be very 
broad resonances. Nevertheless candidates for different continue to exotics exit PQ! Here 
we specialize in tetraquarks, the less difficult multiquarks to compute beyond the baryons 
and hybrids. Notice that there are many possible sorts of tetraquarks: 

- the borromean 3-hadron molecule 

- the Heavy-Heavy-antilight-antilight 

- the hybrid-like tetraquark 

- the Jaffe-Wilczek diquark-antidiquark with a generalized Fermat string 

The borromean 3-hadron molecule 

In an exotic channel, quark exchange leads to repulsion, while quark-antiquark annihila- 
tion is necessary for attraction. A possible way out is adding another meson, allowing for 
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Figure 1: (left) In a hybrid, flux tubes divide into two fundamental flux tunes, one 
connecting the octet with the quark and another connecting the octet to the antiquark. In 
the baryon and in the three-gluon glueball, static quenched Lattice QCD simulations also 
show confinement via fundamental flux tubes, (right) Triple flip-flop Potential potential. 
To the list of potentials to minimize including usually only two different meson pair 
potentials, we join another potential, the tetraquark potential. 



annihilation, to bind the three body system. This has already led to the computation of 
decay widths, which turned out to be wide (21 [3]. 

The Heavy-Heavy-antilight-antilight 

The heavy quarks are easy to bind since the kinetic energy p 2 / (2m) is smaller, thus their 
Coulomb short distance potential could perhaps provide sufficient binding, while the light 
antiquarks would form a cloud around them [I]. 

The hybrid-like tetraquark 

Possibly a quark and antiquark may be in a colour octet, and then the tetraquark is equiv- 
alent to a quark-gluon-antiquark hybrid. Recently we computed in Lattice QCD the color 
fields for the static hybrid quark-gluon-antiquark system, and studied microscopically the 
Casimir scaling [5]. 

Notice that our lattice simulation shows that flux tubes prefer to divide into funda- 
mental flux tubes, or flux tubes carrying a colour triplet flux, as we show in Fig. [TJ 

The Jaffe-Wilczek diquark-antidiquark with a generalized Fermat string 

Since there is no evidence for long distance polarization forces, or Van der Waals forces, 
in hadron-hadron interactions, the two-body confinement potentials cannot be right for 
multiquarks [6]! A solution to this problem consists in considering the flip-flop potential, 
where confining flux tubes or strings take the geometry minimizing the energy of the 
system. Quark Confinement And Hadronic Interactions [TJ. 

Again the flux tubes in the tetraquark are expected to divide and link into fundamental 
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Figure 2: (left) Cherry in a broken glass. Our simplified two-variable toy model is anal- 
ogous to the classical mechanics textbook problem of a cherry in a glass, but a broken 
one where the cherry may escape from. Here we solve this model in quantum mechanics, 
addressing the decay widths of a system compact in one variable and open in the other, 
(right) Plot of our simplified flip-flop potential, as a function of the two radial variables r 
(compact) and p (open). 



flux tubes, and a possible configuration is in a H-like or butterfly-like flux tube. This 
tetraquark can be classified as a Jaffe-Wilczek one since the quarks are combined in a 
diquark-like antitriplet and the antiquarks are combined in a antidiquark-like triplet [8]. 

The technical difficulty in that framework is to compute the decay widths since this 
tetraquark is open for the decay into a pair of mesons. Moreover it is expected that 
the absence of a potential barrier above threshold may again produce a very large decay 
width to any open channel, although Marek and Lipkin suggested that multiquarks with 
angular excitations may gain a centrifugal barrier, leading to narrower decay widths [9]. 

Here we continue a previous work, where we assumed confined (harmonic oscilator- 
like) wavefunctions for the confined objects, one tetraquark and two different pairs of final 
mesons, and computed their hamiltonian. We utilized the Resonating Group Method and 
were surprised by finding very small decay widths [TU] . 

Our approach to study the tetraquark with a generalized Fermat string 

We thus return to basics and decide to have no overlaps. We want to solve the Schrodinger 
equation for the four particles, and from the Schrodinger solutions also compute the decay 
widths. Our starting point is the extended triple flip-flop potential [llj, obtained mini- 
mizing the three lengths depicted in Fig. [TJ Recently, we devised a numerical algorithm 
to compute the Fermat points of the tetraquark and the tetraquark potential \V2\ . 

Solving the Schrodinger equation is then a well defined problem which should be 
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Figure 3: (left) Semi-localized state, or resonance for l r = 1. (right) Bound state for l r = 3. 



solvable, placing our system in a large 12 dimensional box. However this is a very difficult 
problem. Even assuming s-vaves, we would get 3 variables, some confined and some in 
the continuum (similar to problems in extra compactified dimensions or to lattice QCD) 
so we decide to work in a toy model, where the number of variables is simplified. We 
thus simplify the triple flipflop potential, with a single inter-meson variable, using the 
approximation on the diquark and anti-diquark Jacobi coordinates, 



of having a single internal variable p in the mesons. We get a flipflop potential where p 
is open to continuum and r is confined, minimizing only two potentials, 



Our problem is similar to the classical student's problem of a Cherry in a glass. However 
this is not a simple student's problem since the glass is broken and the cherry may escape 
from the glass! The flip-flop and broken glass potentials are depicted in Fig. [2J Here we 
report on our answer [T3] to the question, in the quantum case, are there resonances, and 
what is their decay width? 

2 Finite difference method 

Since there is a single scale in the potential and a single scale in the kinetic energy, we 
can rescale the energy and the coordinates, to get a dimensionless equation, 



Pl3 = P24 





(2) 
(3) 



H${r, p) = [-A r /2 - A p /2 + min(r + V3p, 2r)]$(r, p) = E$(r, p) 



(4) 
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Figure 4: (left) Momenta of the various components as a function of the energy, (right) 
"Phase shifts" obtained from the finite differences ( by projecting the eigenstates in the 
meson- meson eigenstates ). As can be seen the behaviour is irregular when we have more 
than one channel, this is due to the different contributions of multiple channels, for each 
eigenstate calculated in the finite difference scheme. 



that we first solve with the finite difference method. Thus, our results and figures are 
dimensionless. This case is adequate to study equal mass quarks, where the mesons and 
the tetraquark have no constant energy shifts. For instance that would be ok for the light 
tetraquark and meson-meson system 

uudd {s =2) «-> P + P + , (5) 

or the heavy quark system 

CCCC( S =2) ^ J/ip J/lfj. (6) 

We discretize the space in anisotropic lattices and solve the finite difference Schrodinger 
equation, in up to 6000 x 6000 sparse matrices (equivalent to 40 points in the confined 
direction x 150 points in the radial continuum direction). We first look for localized 
states, selecting among the 6000 eigenvalues the ones more concentrated close to the 
origin at p = 0. 

To measure the momenta hi and the phase shifts 5,, we simply fit the large p region of 
the non- vanishing where i indexes the factorized Airy wavefunction in r, the expression 

^ A MhP±3_ (7) 
p 

As can be seen in Fig. HI the momenta ki obey the relation 



h{E) = ^2{E-e i ) ) (8) 

where q is the threshold energy of the respective channel. 

However, the phase shifts we get are not only discrete but rather irregular above 
threshold. In the next Section 3 we compute the phase shifts with an improved method. 
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Figure 5: (left) S-wave scattering cross sections from the channel with l r = and n r = 0. 
(right) S-wave scattering cross sections from the channel with l r = 1 and n r = 0. 



3 Outgoing spherical wave method 

Because the finite difference method is not entirely satisfactory for the computation of 
the phase shifts 5, we move to another method, consisting in in studying the outgoing 
spherical waves. Since the finite difference method shows clearly bands for the different 
internal energies of the mesons, we integrate the confined coordinate r with eigenvalues 
of the meson equation, i.e. with Airy functions, and thus we are left with a system of 
ordinary differential equations in the coordinate p . 

Projecting onto the p coordinate 

We can reduce our problem in the dimensions p, r to a one-dimensional problem in p but 
with of coupled channels. We just have to expand the two-dimensional wavefunction as 

$(r,p) = ]T^(p)<Mr) , (9) 

i 

where the fa are the eigenfunctions of the r confined hamiltonian. The one-dimensional 
potentials are given by 

V l3 {p) = jdh <t>*{r){V FF {r,p)-V MM {r))fa{r) (10) 

where we subtract Vmm from the potential, since Hmm is already accounted for in its 
eigenvalues and eigenfunctions, used for instance in Eq. (j9]). 

Phase shifts 

We now compute the phase shifts, in order to search for resonances in our simplified 
flip-flop model. Solving the outgoing spherical Eq. for this system we can compute the 
partial cross sections and the total cross section for the partial wave I — either directly 
or by using the optical theorem — and determine the phase shifts as well. 
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Figure 6: Comparison of the phase shifts for l r = 0, 1, 2 and 3, with n r = 0. 



Note that our flip-flop potential has the same scales of the simple Schrodinger equation 
for a linear potential, which has a single dimension 



En 



m 



the only energy scale we can construct with h, a and m, the three relevant constants in 
the non-relativistic region. Thus the number of non-relativistic boundstates or resonances 
is independent both of the quark mass m and of the string constant a. 



The centrifugal barrier effect 

Note that we have two distinct angular momenta, which are both conserved, L r = r x p r 
and L p = p x p p . So, each assimptotic state is indexed by its angular momentum l T 
and its radial number n r , and the scattering partial waves are indexed by l p . Thus the 
system can be diagonalized not only in the scattering angular momenta L p but also on the 
confined angular momenta L r . We can describe the scattering process with four quantum 
numbers: The scattering angular momentum l p , the confined angular momentum l r and 
the initial and final states radial number in the confined coordinate r, n.j and rij. 

On Fig. [S]we show the l p — partial cross sections for the scattering from the channel 
with l r = and for l r — 1, with n r = 0. Interestingly, the bumps in the cross section seem 
to occur prior to the opening of a new channel. 

In Fig. Owe compare the phase shifts for different values of l r , namely for l r = 0, 1, 2 
and 3. For l r = 0, we don't observe a resonance, since the phase shift doesn't even cross 
7r/2. However, for the l r = 1 and l r = 2 cases, the phase shifts clearly cross the 7r/2 
line, and a resonance is formed. This behaviour is somewhat expected, since a centrifugal 
barrier in r would, in the case of a true tetraquark, maintain the two diquarks separated, 
favouring the formation of a bound state. The tendency of greater stability for greater 
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6.116 


0.037 
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6.855 


0.131 
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7.462 


0.352 



Table 1: Decay widths function of L. 



orbital angular momenta seems to be further confirmed by the l r — 3, where besides 
the resonance, a true bound state seems to be formed, as can be seen by the different 
qualitative behaviour of the phase shifts for this case. This bound state formation confirms 
our observation of a localized states in Section 2, with the finite difference simulation. 

Finally we can compute the decay width utilizing the phase shift derivative, T/2 = 
(dd/dE)^ 1 computed when the phase shift 5 crosses ir/2, and get the results of Table [TJ 
For instance, for light quarks where m ~ v^I ") — ^00 MeV this results in a l r = 1 decay 
width close to 15 MeV. 

4 Conclusion and outlook to tetraquarks 

We study pentaquarks in the Jaffe-Wilczek model, with a H/butterfly string, but include 
the open channels of decays to meson-meson pairs. We consider an extended flip-flop 
model, where we add the tetraquark string to the two-meson strings. We first apply the 
RGM method assuming that the mesons have gaussian wavef unctions, and we obtain very 
narrow widths. 

We then utilize an approximate toy-model, simplifying the number of Jacobi variables. 
The model is similar to the model of a Cherry in a Broken Glass. This allows the solution 
of the Schrodinger equation with finite differences in a box, where we look for localised 
states, and try to compute phase shifts. 

To compute clearly the phase shifts we then solve the Schrodinger equation for the 
outgoing spherical waves. We compute de decay widths from the phase shifts, and we find 
relatively narrow decay widths. When the produced mesons are unstable, the total decay 
width of the tetraquark is then dominated by the final decays of the produced mesons. 
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